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Abstract — It  has  previously  been  shown  that  dispersion  compen¬ 
sation  can  dramatically  affect  the  interchannel  crosstalk  in  a  mul¬ 
tichannel  analog  fiber  transmission  link.  In  this  work,  we  use  a 
genetic  algorithm  to  find  five-segment  dispersion  maps  that  yield 
low  crosstalk  levels  over  two  octaves  of  microwave-frequency  band¬ 
width  when  amplitude  modulation  is  used.  The  genetic  algorithm 
suggests  that  optimal  dispersion  maps  have  low  residual  disper¬ 
sion.  Despite  the  genetic  algorithm’s  ability  to  optimize  dispersion 
maps  with  many  fibers,  it  is  possible  to  obtain  similar  crosstalk 
levels  from  a  simpler  two-segment  design  whose  dispersion  is  fully 
compensated. 

Index  Terms — Analog  transmission,  cross-phase  modulation 
(XPM),  crosstalk,  intensity  modulation,  self-phase  modulation, 
wavelength-division  multiplexing. 


I.  Introduction 

IN  wavelength-division-multiplexed  fiber-optic  analog  sys¬ 
tems,  cross-phase  modulation  (XPM)  causes  amplitude 
modulation  (AM)  in  one  channel  to  be  transferred  to  phase 
modulation  (PM)  in  the  others,  and  chromatic  dispersion  then 
converts  this  PM  into  undesired  AM  [1],  [2].  Unlike  digital 
systems,  where  the  input  optical  powers  are  low  and  the  signal 
is  periodically  reamplified,  somewhat  reducing  this  effect, 
analog  systems  often  employ  high  input  optical  powers  with  no 
reamplification  to  obtain  a  high,  spurious-free  dynamic  range 
[3]  over  a  bandwidth  of  several  octaves.  As  a  result,  the  effect 
of  interchannel  nonlinear  crosstalk  can  be  much  more  severe 
in  these  analog  systems  [2],  [4].  Several  authors  have  pre¬ 
sented  work  investigating  the  impact  of  interchannel  crosstalk, 
and  several  have  presented  novel  methods  for  reducing  it  by 
modifying  the  modulation,  detection,  or  transmission  scheme 
[5]— [7] .  However,  these  methods  either  require  complicated 
and  costly  receivers,  require  parallel  fiber  links,  or  work  only 
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over  a  limited  bandwidth.  Additionally,  simple  methods  for 
reducing  the  crosstalk,  such  as  polarization  multiplexing,  offer 
only  modest  reductions  in  XPM. 

The  role  of  chromatic  dispersion  in  converting  PM  due 
to  XPM  into  intensity  modulation  suggests  that  tailoring  the 
dispersion  map  of  a  link  or  delay  line  by  splicing  together 
several  fibers  with  different  dispersion  values  may  be  a  useful 
and  simple  technique  to  reduce  crosstalk  in  analog  systems. 
Although  some  authors  have  investigated  using  periodic  dis¬ 
persion  compensation  to  reduce  XPM  in  long  or  installed  links 
[8]— [10],  little  work  has  been  presented  investigating  the  possi¬ 
bility  of  designing  a  fiber  span  yielding  low  crosstalk  using  an 
arbitrary  dispersion  map.  In  this  letter,  we  use  both  a  genetic 
algorithm  for  a  five-segment  design  and  a  global  search  for  a 
two-segment  design  to  explore  dispersion  maps  that  lead  to  low 
interchannel  crosstalk  over  a  range  of  modulation  frequencies. 
We  show  that  while  the  choice  of  the  fiber  segment  lengths  and 
dispersions  can  significantly  impact  the  maximum  crosstalk 
level,  the  best  five-segment  system  offers  no  significant  advan¬ 
tage  over  an  optimized  two-segment  system. 

II.  Methods 

Previous  work  shows  that  dispersion  compensation  can  be 
used  to  reduce  crosstalk  at  individual  modulation  frequencies 
[4] .  Here,  we  simulated  systems  comprising  concatenated  fiber 
segments,  each  with  constant  dispersion,  and  we  used  a  ge¬ 
netic  algorithm  to  find  configurations  that  lead  to  small  residual 
crosstalk  over  a  desired  range  of  modulation  frequency.  The  ge¬ 
netic  algorithm  attempts  to  find  a  global  optimum,  so  that  we 
have  some  assurance  that  the  crosstalk  value  corresponding  to 
the  configuration  found  by  the  algorithm  is  close  to  the  best 
achievable  value. 

We  initiate  the  genetic  algorithm  simulation  by  creating 
a  randomly  chosen  population  using  the  constraints  of  the 
problem — a  process  known  as  genesis.  In  our  simulation  of 
the  crosstalk,  each  individual  is  a  data  structure  consisting 
of  the  lengths  and  local  chromatic  dispersion  values  of  the 
concatenated  fiber  segments.  We  then  evaluate  the  fitness  of 
each  individual  by  computing  the  negative  of  the  maximum 
crosstalk  in  the  target  range  of  modulation  frequency.  We  use 
the  negative  of  the  maximum  crosstalk  so  that  lower  crosstalk 
corresponds  to  higher  fitness.  Once  the  fitness  of  each  member 
of  the  population  is  evaluated,  we  rank  the  members  and  use  a 
selection  process  to  create  the  next  generation  of  the  popula¬ 
tion,  favoring  individuals  with  high  fitness  over  those  with  low 
fitness.  The  information  in  the  data  structure  describing  each 
individual  is  then  exchanged  among  randomly  chosen  pairs  of 
individuals — crossover — and  then  a  mutation  process  ensures 
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that  new  randomly  chosen  parameter  values  are  introduced 
into  the  population.  After  the  crossover  and  mutation  steps  are 
complete,  the  next  generation  is  evaluated  for  fitness,  a  new 
selection  process  is  performed,  and  the  process  repeats  until  the 
crosstalk  level  is  sufficiently  improved. 

In  the  genesis  procedure,  we  choose  a  population  of  100  indi¬ 
viduals,  whose  system  parameters  are  restricted  in  the  following 
ways.  The  segments’  chromatic  dispersion  values  can  take  any 
value  in  the  range  of  1  to  50  ps/nm-km  with  either  positive  or 
negative  sign.  The  total  length  of  the  system  is  kept  fixed  at 
15  km,  an  intermediate  length  for  a  wideband  analog  system, 
while  the  lengths  of  each  of  the  fiber  segments  are  allowed  to 
be  any  fraction  of  the  total  length.  In  our  simulations,  we  fix 
the  fiber  loss  at  0.2  dB/km  and  the  nonlinear  coefficient  7  at 
1.55  W-1  •  km-1  for  all  segments  in  the  system.  We  also  use 
an  average  launch  power  of  10  mW  per  channel  (the  launch 
power  required  to  produce  a  spur-free  dynamic  range  of  ap¬ 
proximately  110  dB-Hz2/3  for  this  link)  and  the  same  200-GHz 
channel  spacing  used  in  our  previous  work  [1],  [4]. 

Next,  we  evaluate  the  fitness  of  each  individual  by  inte¬ 
grating  the  system  of  four  ordinary  differential  equations 
(ODEs)  derived  in  [1]  for  values  of  the  modulation  frequency 
in  the  two-octave  range  of  4-16  GHz  with  a  1-GHz  step  size. 
We  use  this  method  because  solving  the  ODE  system  is  signif¬ 
icantly  computationally  faster  than  solving  the  full  nonlinear 
Schrodinger  partial  differential  equation.  After  evaluating  the 
fitness  of  each  individual,  the  population  is  sorted  by  fitness 
value. 

Once  the  population  is  sorted,  we  use  a  roulette  wheel  se¬ 
lection  process  to  determine  the  next  generation  [11].  In  this 
process,  each  individual  in  the  population  is  assigned  a  wedge 
on  a  virtual  roulette  wheel.  The  angle  subtended  by  each  wedge 
is  proportional  to  the  individual’s  fitness  value  minus  the  min¬ 
imum  fitness  of  the  population,  so  that  the  least  fit  individual 
is  eliminated.  A  set  of  100  independent  random  numbers  uni¬ 
formly  distributed  between  0  and  27t  are  then  chosen,  corre¬ 
sponding  to  spins  of  the  roulette  wheel,  and  the  wedge  corre¬ 
sponding  to  each  spin  selects  the  corresponding  individual  to 
survive  for  the  next  generation.  In  this  process,  it  is  likely  to 
have  multiple  copies  of  the  most  fit  individuals. 

After  selection,  we  perform  crossover,  in  which  randomly  se¬ 
lected  pairs  of  individuals  exchange  parameters  of  their  fiber 
segments.  After  pairing,  we  choose  a  random  number  between 
zero  and  one  and  perform  crossover  on  the  pair  if  the  number  is 
less  than  0.5.  If  crossover  occurs,  we  randomly  select  a  fiber  seg¬ 
ment,  referred  to  as  the  crossover  point,  within  one  of  the  indi¬ 
viduals.  The  lengths  and  dispersion  values  of  the  fiber  segments 
beyond  the  crossover  point  are  exchanged  between  the  pair.  We 
have  chosen  to  constrain  the  system  length  to  remain  constant 
at  15  km.  However,  after  the  individuals  exchange  their  param¬ 
eters,  the  system  length  will  usually  be  longer  or  shorter  than 
15  km.  Therefore,  after  crossover,  we  change  the  lengths  of  the 
segments  by  simply  rescaling  the  system  back  to  15  km.  In  per¬ 
forming  this  rescaling,  we  change  the  local  dispersion  values  so 
that  the  accumulated  dispersion  in  each  fiber  segment  is  main¬ 
tained.  In  this  way,  the  crossover  process  preserves  the  accumu¬ 
lated  dispersion  of  all  fiber  segments  in  the  population,  although 
their  lengths  are  changed. 


Distance  (km)  Modulation  Frequency  (GHz) 


Fig.  1.  (a)-(c)  Five-fiber  dispersion  maps  resulting  from  the  top  three  popu¬ 
lation  members  after  15  generations  of  the  genetic  algorithm,  and  (d)-(f)  the 
corresponding  crosstalk  in  decibels  as  a  function  of  modulation  frequency.  The 
fitness  (negative  maximum  crosstalk  in  4-16  GHz)  increases  from  (a)  to  (c). 

Finally,  in  the  mutation  process,  we  select  a  random  number 
between  zero  and  one  for  each  fiber  segment  in  the  population. 
If  the  number  is  less  than  0.15,  we  perform  a  randomization  of 
the  local  dispersion  value  for  that  fiber  segment,  while  keeping 
its  length  constant.  The  new  dispersion  value  is  again  chosen 
to  be  in  the  range  of  1-50  ps/nm-km  with  either  a  positive  or 
negative  sign.  While  the  crossover  process  is  devised  to  main¬ 
tain  the  accumulated  dispersion  while  changing  the  segments’ 
lengths,  the  mutation  process  maintains  the  segments’  lengths 
while  changing  the  accumulated  dispersion. 

After  several  generations,  the  genetic  algorithm  eliminates 
configurations  with  high  maximum  crosstalk  in  favor  of  those 
with  low  crosstalk.  Our  algorithm  finds  examples  of  configura¬ 
tions  with  very  low  crosstalk  levels  over  the  range  of  4-16  GHz. 

We  computed  15  generations  of  the  population  using  the  ge¬ 
netic  algorithm.  We  verified  that  the  algorithm  works  by  noting 
that  the  ensemble  average  fitness  level  increases  as  a  function  of 
the  generation  number.  After  15  generations,  the  average  fitness 
level  does  not  increase  appreciably  between  successive  gener¬ 
ations.  In  Fig.  1,  we  show  the  dispersion  map  and  the  corre¬ 
sponding  crosstalk  curves  for  the  top  three  individuals.  The  dis¬ 
persion  map  parameters  shown  in  Fig.  l(a)-(c)  look  very  dif¬ 
ferent  from  each  other. 

Despite  these  differences,  the  maximum  crosstalk  achieved  in 
each  of  the  three  cases  is  close  to  —40  dB.  The  first  generation 
of  the  genetic  algorithm  yielded  a  similar  result  for  the  lowest 
maximum  crosstalk  level.  The  genetic  algorithm  was  only  mar¬ 
ginally  successful  at  lowering  this  best  crosstalk  level  in  future 
generations. 

Several  salient  points  arise  from  these  data.  First,  the  disper¬ 
sion  map  can  indeed  be  tailored  to  lower  the  crosstalk  level  over 
a  desired  range  of  frequencies.  Second,  finding  an  example  of 
a  “good”  dispersion  map  is  not  difficult  and  does  not  require  a 
complicated  optimization  procedure.  Finally,  while  tailoring  the 
dispersion  map  is  beneficial,  lowering  the  crosstalk  level  fur¬ 
ther  will  require  changes  in  the  transmitter  or  receiver  design. 
In  sum,  it  is  unlikely  that  there  is  a  special  dispersion  map  that 
is  dramatically  better  than  all  the  others.  Many  dispersion  maps 
produce  nearly  the  same  improvement. 

Fig.  2  shows  a  scatter  plot  of  the  maximum  crosstalk  for  all 
100  individuals  in  each  of  the  last  three  generations  of  the  ge¬ 
netic  algorithm  as  a  function  of  the  residual  dispersion  in  the 
line.  The  most  fit  individuals  appear  to  be  those  that  have  a 
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Residual  dispersion  (ps/nm) 


Fig.  2.  Scatter  plot  showing  the  relationship  between  the  maximum  crosstalk 
and  the  residual  dispersion  for  all  population  members  in  the  last  three  (13-15) 
generations  of  the  genetic  algorithm.  Solid  circles  correspond  to  generation  13, 
diamonds  correspond  to  generation  14,  and  squares  correspond  to  generation  15. 


Fig.  3.  Contour  plots  of  the  maximum  crosstalk  in  the  4-  to  16-GHz  range 
as  a  function  of  compensator  length  and  fiber  dispersion  D  for  a  100% 
compensated  two-fiber  dispersion  map.  The  total  length  of  the  dispersion  map 
(fiber  +  compensator)  was  15  km.  In  (a),  the  compensator  comes  before  the 
fiber  (pre-compensation),  and  in  (b)  the  compensator  comes  after  the  fiber 
(post-compensation). 


crosstalk  correspond  to  configurations  with  low  dispersion  in 
the  link  fiber  and  a  compensator  length  that  is  the  same  as  the 
link  fiber.  Finally,  a  post-compensated  link  outperforms  a  pre¬ 
compensated  link.  This  difference  is  due  to  the  interaction  of 
the  fiber  dispersion  and  nonlinearity  in  the  link.  In  all  cases,  the 
maximum  power  in  the  link  is  at  the  launch  point,  and  the  sign 
of  the  dispersion  of  the  first  fiber  determines  the  interaction  of 
the  dispersion  and  the  nonlinearity.  Currently,  the  reason  for  the 
better  performance  of  the  post-compensated  configuration  is  not 
known. 


IV.  Conclusion 

We  used  a  genetic  algorithm  to  find  five-segment  dispersion 
maps  that  minimize  the  crosstalk  over  two  octaves  of  bandwidth 
when  AM  is  used  in  an  analog  fiber-optic  link.  We  found  that 
a  large  number  of  dispersion  maps  with  widely  varying  map 
parameters  yielded  nearly  the  same  advantage — at  best  about 
—40  dB  for  our  15-km  system.  We  also  found  that  it  is  possible 
to  achieve  almost  the  same  advantage  using  two-segment  maps 
with  a  zero  net  dispersion. 

Our  conclusion  is  that  a  significant  advantage  can  be  obtained 
by  using  appropriate  dispersion  maps.  However,  there  is  little 
or  no  advantage  to  using  complex  dispersion  maps.  It  does  not 
appear  possible  to  significantly  outperform  a  properly  chosen 
two-segment  dispersion  map. 


residual  dispersion  that  is  close  to  zero.  Therefore,  it  is  natural 
to  restrict  our  search  for  optimal  dispersion  maps  to  those  that 
have  zero  net  dispersion. 

III.  Simulations  of  Two-Segment  Dispersion  Maps 

Motivated  by  the  results  of  Section  II,  we  next  investigated 
dispersion  maps  comprised  of  two  segments — a  link  segment 
and  a  compensating  segment.  As  before,  the  total  system  length 
is  restricted  to  15  km.  The  compensating  segment  compensates 
fully  for  the  dispersion  in  the  link  segment  so  that  the  residual 
dispersion  is  zero.  We  considered  two  cases — a  pre-compen- 
sated  system  and  a  post-compensated  system.  In  each  case,  the 
length  of  the  compensator  is  varied  between  0.5  and  7.5  km,  and 
the  link  fiber  dispersion  is  varied  between  6  and  27  ps/nm-km. 
The  attenuation  in  each  segment  was  set  to  0.2  dB/km. 

Fig.  3  shows  contour  plots  of  the  maximum  crosstalk  value 
achieved  within  the  two-octave  range  of  modulation  frequencies 
of  4-16  GHz.  In  Fig.  3(a),  the  system  is  pre-compensated,  while 
in  Fig.  3(b),  the  system  is  post-compensated.  In  our  simulations, 
we  integrate  the  ODE  system  described  above  for  the  range  of 
4-16  GHz  with  a  1-GHz  resolution  and  interpolate  quadratically 
to  find  the  maximum  crosstalk. 

Fig.  3  shows  that  a  two-segment  dispersion  map  can  achieve  a 
similar  crosstalk  level  as  the  five-segment  dispersion  maps  that 
were  investigated  using  the  genetic  algorithm  (about  —42  dB 
in  these  simulations).  Additionally,  when  using  either  pre-  or 
post-compensation,  the  configurations  that  generate  the  smallest 
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